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Abstract — Protocol sequences are used in channel access for 
the multiple-access collision channel without feedback. A new 
construction of protocol sequences with a guarantee of worst- 
case system throughput is proposed. The construction is based 
on Chinese remainder theorem. The Hamming crosscorrelation 
is proved to be concentrated around the mean. The sequence 
period is much shorter than existing protocol sequences with the 
same throughput performance. The new construction reduces the 
complexity in implementation and also shortens the waiting time 
until a packet can be sent successfully. 

Tags: Protocol sequences, collision channel without feed- 
back, wobbling sequences. 

I. Introduction 

D Protocol sequences are periodic binary sequences for 
multiple-access control in the collision channel without feed- 
back 1 6 1. In a time-slotted scenario, each user repeatedly reads 
out the value of a statically assigned protocol sequence, and 
sends a packet in a time slot if the sequence value is equal 
to one. If two or more users transmit simultaneously in the 
same time slot, there is a collision and the collided packets 
are assumed unrecoverable. If there is exactly one transmitting 
user while the others remain silent, the received packet is 
assumed error-free. We do not assume any feedback from 
the receiver and any coordination among the transmitters. 
This assumption is applicable to low-cost and low-complexity 
wireless sensor networks, as it is not necessary to spare any 
hardware on monitoring the channel; the transmitters simply 
send a packet whenever the value of the assigned protocol 
sequence is one, regardless of the channel condition. For 
simplicity in presentation, we assume slot-synchronization in 
this paper This requirement can be relaxed without much 
degradation in performance. 

This channel model is considered in practical sensor net- 
works, such as f-MAC 18]. One of the main design issues 
is the construction of protocol sequences. Some design ob- 
jectives are addressed in |9|. Since we do not assume any 
coordination among the users, the users may not start their 
protocol sequences at the same time. This incurs relative delay 
offsets among the transmitters. Our first objective is to design 
protocol sequences such that no matter what the relative delay 
offsets are, the system throughput is provably larger than some 
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positive constant. This provides a throughput guarantee in the 
worst-case. 

The second objective is to minimize the sequence period, 
which measures the delay one has to wait until the promised 
number of successful packets go through the channel. For 
protocol sequences with very long period, a user may suffer 
starvation in the short term even though the throughput in 
the whole period is very good. This issue is alleviated if the 
sequence period is minimized. 

The two objectives mentioned above are contradicting, and 
there is a tradeoff between them. In Q, |6l, f9l, a class of pro- 
tocol sequences, called shift-invariant sequences are studied. 
This class of protocol sequences can achieve optimal system 
throughput, but the sequence period grows exponentially as a 
function of the users, and hence is not of practical interests 
when the number of users is large. Another class of protocol 
sequences, called wobbling sequences is constructed in ifTOl . 
The system throughput is provably larger than 0.25 for any 
choice of relative delay offsets, and the sequence period grows 
like M^, where M is for the number of users in the system. 
In this paper, we construct shorter protocol sequences with 
roughly the same throughput performance by the wobbling 
sequences. 

Other constructions of protocol sequences are investigated 
in m, Q, sometime under the name of cyclically permutable 
constant weight codes (CPCWC). The difference between 
CPCWC and our protocol sequences is that the latter only 
has Hamming crosscorrelation requirement, but the former has 
both autocorrelation and crosscorrelation constraints. In 1*21, 
constructions using optical orthogonal codes and cyclic super- 
imposed codes are considered. 

In Section [III the new construction of protocol sequences 
is described. In Section [Till we investigate the crosscorrelation 
properties, which are crucial in the derivation of a lower bound 
on system throughput in Section HV] Comparison with shift- 
invariant and wobbling sequences is given in Section ITV] 

II. Construction of Protocol Sequences via 
Chinese Remainder Theorem 

We will use sequence "period" and sequence "length" 
interchangeably. Let Z„ = {0, 1, . . . , n — 1} be the residues 
of integers mod n. The components in a sequence of length L 
are indexed by Z/,. The Hamming weight of a binary sequence 
a{t) of length L is the number of ones in a{t) in a period. For 



two binary sequences a{t) and b{t) of length L, their Hamming 
crosscorrelation function is defined as 

Hab{T) := <^it)Kt + r), 

where t is the delay offset. 

We shall construct sequences with length L = pq, where 
p and q are two relatively prime integers. In this paper, we 
will take p to be a prime number and q an positive integer not 
divisible by p. By the Chinese Remainder Theorem (CRT) ijS), 
there is a bijection between Zp^ and the direct sum 

Gp^q '.— Zp © • 

The bijective mapping $ : Zpg -> Gp^g is given by 

$(a;) — {x mod p, x mod q). 

Henceforth, we will identify Zp^ with Gp.g. 

CRT Construction Given a prime number p and an integer 
q relatively prime to p, we define a sequence Sg it) of length 
pq, for g = 0,1,2. . 1, by 

"5(0 = 0'5,j) for some j, < j < q, 

Sg(t) := < 

I otherwise, 

where j is the residue of j in Zp. We call the sequence Sg{t) 
the CRT sequence generated by g. The integer g is called the 
generator of Sg(t). 

Alternately, we can define the CRT sequences by specifying 
their characteristic sets. For each g G {0, 1, . . . ,p — 1}, let 

Tg:={{]g,j)eGp^q: 0<j <q}. (1) 

We note that Ig is an arithmetic progression in Gp^q with 
common difference {g, 1). The CRT sequence with generator 
g is obtained by setting Sg{t) = 1 if and only if G Ig. 
The Hamming weight of each sequence is equal to q. 

Example 1. p = 5 and q = 9. The five CRT sequences so{t) 
to S4{t) are listed as follows. The common period is 45, and 
the Hamming weight of each sequence is equal to 9. 

50 : 10000 10000 10000 10000 10000 10000 10000 10000 10000, 

51 : 11111 11110 00000 00000 00000 00000 00000 00000 00000, 

52 : 10000 10000 00010 00000 01000 01001 00001 00100 00100, 

53 : 10000 10000 01000 01000 00100 00010 00010 00001 00001, 

54 : 10000 10000 00100 00101 00001 00000 01000 00010 00010. 

III. Crosscorrelation Properties 

The main idea of using CRT in constructing protocol 
sequences hinges on the fact that $ is a homomorphism of 
abelian groups, so that the analysis of crosscorrelation can be 
carried out in Gp^q instead of Zp^. We remark that in fQ, 0, 
Q, HI, the idea of CRT appears in the same way as in this 
paper. 

Given a one-dimensional delay r, we denote its two- 
dimensional counterpart by (ti,T2) := $(t). For h G Zp, 
define the translation Ih + (ti, T2) of Ih by (ti, T2) by 

{{x + Ti,y + T2) e Gp^q : {x,y)£lh}. 



The Hamming crosscorrelation Hgh{T) can be computed by 

where |^| indicates the cardinality of a set A. To distinguish 
arithmetic operations in Zp and Z^, we use ©p and 0p for 
addition and subtraction in Zp, and (Bq and Qq for addition 
and subtraction in Zq. In this notation, we have 

Ih + (ri, T2) = {{jh ®p n, j ®g r2) : j = 0, 1, . . . , q - 1}. 

By a change of variable, Ih + (ti,T2) can be written as 

{(((jegT2)/i) ®pTi,j) : j = 0,l,...,(Z-l}. 

After comparing with the definition of Ig in ([U, we see that 
\Ig n {Ih + {ti,T2))\ is equal to the number of solutions to 

xg = {{x Qq T2)h) ©p Ti mod p. (2) 

for X = 0,1, . . . , q — 1. The problem of computing the cross- 
correlation function is thus reduced to counting the solutions 
to ©. 

The following simple lemma is useful in the derivation of 
Hamming crosscorrelation. 

Lemma 1. Let p be a prime number For each b G Zp, the 
number of solutions to x = h mod p for x going through d 
consecutive integers c, c+l,...c + d— 1, equals 

d/p if p divides d, 

\d/p\ + 5 otherwise, 

where 5 is either or 1. 

Proof: In the first case where d is divisible by p, if we 
reduce the integers c, c+\, . . . ,c + d—l mod p, we have each 
element in Zp repeated exactly d/p times. Hence, for each 
b G Zp, there are exactly d/p integers in {c, . . . , c + d — 1} 
whose residue mod p equal 6 

For the second case, where d is not divisible by p, we divide 
the d consecutive integers into two parts. The first part consists 
of \d/p\p consecutive integers and the second part consists of 
the last d — Yd/p\p integers. Among the first \d/p\p integers, 
exactly [d/pj of them equal b mod p. The residues of the 
d — \ d/p\p integers in the second part are distinct, and hence 
at most one of them is equal to b. The number of integers in 
{c, . . . , c + d — 1} whose residue equal b is either [d/pj or 
\_d/p\ +1. ■ 

From Lemma [T] we obtain the crosscorrelation between 
so(0 ^iid other sequences. 

Theorem 1. For g ^ 0, the Hamming crosscorrelation of 
Sg(t) and so(t) is equal to either \_q/p\ or [g/pj + 1. 

Proof: If we put /i = in dU, we get x = g~^Ti mod 
p. The number of integers in {0,l,...,q — 1} which equal 
g^^Ti modp is either [q/p\ or [q/p\ +1 by Lemma [T] ■ 
For nonzero h, we can divide both sides of (|2|i by /i and 
re-write it as 

x{h^^g) = {x Qq T2) ©p {h^^Ti) modp. 



For each fixed T2, as ti runs through Zp, h^^Ti also runs 
through the complete set of residues mod p. Therefore, the 
distribution of Hamming crosscorrelation between Sg{t) and 
Sh{t) is the same as the distribution of Hamming cross- 
correlation between Sg/h{t) and si{t). We henceforth focus 
on the case h = 1 without any loss of generality. 

To aid the derivation of the Hamming crosscorrelation, we 
first prove the following lemma. 

Lemma 2. Let g € Zp\{l}, and denote $(t) by (ti,T2). The 
Hamming crosscorrelation between Sg(t) and Si{t), namely 
Hgi^r), satisfies the following properties: 

1) Hgi{T) equals the number of solutions to 

X = ag{Ti,T2) + hgl{Q < X < T2) raodp, (3) 

for X = 0, 1, . . . , g — 1, where 

ag{n,T2):^{g~ir\n~f2), (4) 
hg := {g - l)-^q, (5) 

and I is the indicator function defined as 

, , 1 if P is true, 
I(P) < 

10 if P is false. 

2) Let T and t' denote two relative delay offsets. Suppose 
that the first component of^{T) and $(t') are the same, 
and the second component of ^(t) and ^{t') defer by 
a multiple of p, then then Hgilr) = iJgi(T'). 

Proof: After setting /i in (|2]i to 1, we obtain 

xg = (x Qg T2) (Bp Ti mod p. (6) 

We want to show that the number of solutions to for q = 
0, 1, . . . , g — 1, is the same as the number of solutions to (O. 

We consider x in two disjoint ranges: (i) < x < T2, and 
(ii) r2 < a; < g. In the first case, x Qg T2 is congruent to 
X + q — T2 mod q. So, for < x < T2, ^ is equivalent to 

xg = X + q — f2 + Ti mod p (7) 

where q and f2 are residues of q and T2 in Zp, respectively. 

In the second case, for x — T2, T2 + 1, . . . , g — 1, © is 
equivalent to 

xg = X — f2 + Ti mod p. (8) 
We combine ^ and ^ in one line as 

x{g — 1) = -f2 + Ti + ql{0 < X < T2) mod p, 

Since g is not equal to 1 by assumption, we can divide by 
(g — 1) and obtain (O. This proves the first part of the lemma. 

The second part of the lemma is vacuous if q < p. So 
we assume q > p. (The case q — p is excluded because it is 
assumed that q is relatively prime with p.) Let (ri, T2) — <&(t) 
and {t[,T2) = It is sufficient to prove the statement for 

Ti = t[ and T2 = T2 + p, namely, the number of solutions 
to (O and the number of solutions to 

X = ag(Ti, Tj) + 6gl(0 < X < T2) mod p (9) 



for a; = 0, 1, . . . , q — 1, are the same. We note that ag(Ti, T2) 
is equal to ag(ri,r2). However, the arguments inside the 
indicator function are different. We divide the range of x into 
three disjoint parts: 

Xi :={0,l,...,r2-l}, 
X2 := {T2,r2 + 1, . . . ,r2 + J3 - 1}, 
:= {t2 + p,T2 + p + 1, . . . ,q - 1}. 

Since q > p, is non-empty. For x Xi, I{0 < x < 
T2) = 1(0 < X < T2). Therefore ^ and (|9]l have the same 
number of solutions for x in Xi. For x ^ X2, both (O and 
(|9|l have exactly one solution by Lemma [T] For x £ X3, (O is 
equivalent to (|9]l, and hence has the same number of solutions 
as ^ does. In conclusion, the number of solutions to (OJ and 
(|9]l for X £ U A:2 U ^3 are the same. This finishes the proof 
of the second part of the lemma. ■ 
From now on, we assume that q > p, which is the case of 
practical interest. 

Theorem 2. Let p and q be positive integers such that p is 
prime, gcd(p, q) = 1 and q > p. Let m denote the quotient of 
q divided by p, i.e., m = lq/p\, and let g G Zp, ^ g =i 1. Let 
q be the residue of q mod p, and bg be defined as in (|5]l. The 
Hamming crosscorrelation between Sg{t) and Si{t) is bounded 
between 

m — 1 and m + 1 if < bg < p — q, or (10) 
m and m + 2 if p — q < bg < p. (11) 

Proof: By the second part of the previous lemma, we only 
need to consider T2 = 0, 1, . . . ,p — 1. In this proof, we will 
denote i/gi(r) by Hgi{Ti,T2), with {ti,T2) equal to ^(t). 

We first prove the first case in (fTol i by considering two cases. 

Case 1, < T2 < g 

Suppose that (O has no solution for < a; < T2. As the 
indicator function in (|3]l is zero for a; = T2, T2 + 1, . . . , q — 1, 
^ is reduced to 

X = ag(Ti, T2) mod p. 

The number of integers in {t2, T2 + 1, . . . , q — 1}, say d, 
satisfies [d/p\ — m. By LemmalT] we have either m or m+1 
solutions to ([3]) for a: > r2. 

Secondly, suppose that ([3]) has exactly one solution for < 
a; < T2. The indicator function in (O is equal to 1 for < 
a; < T2. Hence, 

<ag{Ti,T2) + bg <T2. (12) 

We claim that (O has no solution for x — T2,T2 + 1, ■ ■ . ,q—l- 
Otherwise, we have 

T2 < ag{Tl,T2) < q, 

which, after combining with the assumption that 1 < &g < 
p — q — 1, yields 

T2 < ag{Tl,T2) +bg <p-l. 



This contradicts with ( fT2] i and proves the claim. For 

q<x<q, 

there are exactly m solutions by Lemma [T] The total number 
of solutions to (O for x — 0,1, . . . ,q — l, is thus m + 1. Hence 

Hgl{Tl,T2) = m + 1. 

Case 2: q < T2 < p 

By Lemma [1] (|3]l has either or 1 solution for < x < r2, 
and either to — 1 or to solutions for T2 < x < q. Hence, 
Hgi{Ti,T2) is within the range of {to — 1, m, to + 1}. 

For bg = p — q + 1, . . . ,p — 1, we again consider two cases. 
Case 1: < r2 < g 

By Lemma [1] ^ has either or 1 solution for < x < T2, 
and either to or to + 1 solutions for T2 < x < q. Therefore, 

Hgi{Ti,T2) G {m, ni + 1, m + 2}. 
Case 2: q < T2 < p 

Suppose that (O has no solution for < a; < r2, i.e., 

T2 < ag{Tl,T2) +bg < p. (13) 

We claim that ^ must have one solution for x in the following 
range 



T2 < X < p + q- 



(14) 



From the assumption of q < T2 < p, we deduce that 

q<p + q-T2<p, 

so that the range in (fl4l i is non-empty and consists of no 
more than p integers. If the claim were false, we would have 
no solution to (|3]l for T2 < x < p + q, implying that 



q < ag{Tl,T2) < T2. 



(15) 



Here, we have used the fact that the indicator function in ^ 
is equal to zero for x in the range in (fT4l l. By adding (flST l to 

p-q + l<hg<p-l 

and reducing mod p, we obtain 

1 < ag{Ti,T2) +bg < T2, 

which is a contradiction to ( fT3] l. Thus, the claim is proved. 
For X = p + q,p + q + l,...,q — 1, there are exactly m — 1 
solutions to ^ by Lemma [T] Totally there are m solutions, 
and thus Hgi{Ti,T2) — to. 

Finally suppose that (|3]l has exactly one solution for < 
a; < r2 . As the number of solutions to (O for a; = T2 , T2 + 
1, . . . , q — 1 is either m — 1 or m by Lemma[T] the total number 
of solutions to ([3]) is either m or to + 1. 

In any case, we see that Hgi{Ti, T2) is either to, to + 1 or 
TO + 2. ■ 

Thereom |2] asserts that for any pair of distinct CRT se- 
quences, the Hamming crosscorrelation is either between to— 1 
and TO + 1, or between to and to + 2. For the whole sequence 
set, the Hamming crosscorrelation is therefore four-valued. 
We next show that for some special choice of q, namely 
q = ±1 modp, the Hamming crosscorrelation of the whole 
sequence set assumes only three distinct values. 



Theorem 3. Let p and q be integers as in Theorem |2] 

\) If q is of the form mp + 1 for some positive integer to, 
then for g = 2,3, . . . ,p — 1, Hgi (r) is between to — 1 
and TO, + 1. 

2) If q be of the form mp+ {p— I) for some positive integer 
TO, then for 17 = 2, 3, ... ,p — 1, Hgi (r) is between m 
and TO + 2. 

Proof For he first part of the theorem, we have q equal 
to 1 mod p. So 



bg = {g-l)-'q={g-l)-' mod 



P- 



Since g is between 2 and p — 1 inclusively, g — 1 is between 
1 and p — 2, and hence the inverse of g — 1 mod p is also 
between 1 and p — 2. We thus obtain < bg <p— \=p — q. 
The result now follows from Theorem ID 

The second part can be proved similarly from Theorem |2] 
by putting q^ p — \. ■ 

Together with Theorem [T] which says that the Hamming 
crosscorrelation between so(t) and Sg{t), for g ^ 0, is 
either m or to + 1, we prove that the Hamming crosscor- 
relation of the whole CRT sequence set is three-valued when 
q = ±1 mod p. The sequences in Example 1 are generated 
with q = — 1 mod p. We can verify that the Hamming 
crosscorrelation in Example 1 is either 1, 2 or 3. 

IV. Lower Bound on System Throughput 

The three-valued result in Theorem |3] suggests that the 
variation of Hamming crosscorrelation due to relative delay 
offsets is minimal when q = ±1 mod p. We single out the 
q = —1 modp case below and derive a lower bound on 
the resulting system throughput. The case of g = 1 mod p 
is similar and omitted. 

When q is of the form kp — 1 for some integer k > 2, 
the CRT construction yields p protocol sequences of length 
L = kp^ — p and Hamming weight kp — 1. By Theorem |3] 
the largest Hamming crosscorrelation value is fc + 1. 

We pick AI sequences and form the CRT sequence set of 
size p in order to support M users. Here, M is an integer 
whose value will be optimized later. Since a user sends kp—1 
packets in a period, and each other user may collide with him 
in at most fc + 1 packets, the number of successful packets per 
user per period is no less than kp~l — {M — l)(fc + 1). The 
total number of successful packets, summed over all M users, 
is thus lower bounded by 

M[fcp-l-(A-f-l)(fc + l)]. (16) 

By the method of completing square, we can write ( fTSl l as 

fc(p+l)^2 /fc(p+l) 



(fc+1) 



M 



2(fc + l)/ V2(fc + 1 
We see that the maximum value in ( fTSI ) is obtained when 



(17) 



M* 



fc(p+l)/(2(fc + l)). 



(18) 



Since M must be an integer, after taking the floor of ( fTsT l. 
we obtain the following theorem. 



Theorem 4. Let p be a prime number, k > 2, and M' be the 
largest integer smaller than or equal to M* in ( llSl l. By picking 
M' sequences from the CRT construction with parameters p 
and q ~ kp — 1, the system throughput is lower bounded by 



p{kp - 1) 



(fc + 1) 



(19) ^ 



Proof: Consider the expression in ( fTSI l as a function of 
M, and denote it by f{M). Since the difference between M* 
andM'isatmostone, /(A/*)-/(M') < (fc+l)(A/'-Af*) < 
+ After division by the period p{kp — 1), we have the 
following lower bound on system throughput, 

/(M') ^ f(M*) - (fc + 1) 



p{kp - 1) p{kp - 1) 

which can be readily seen to be the same as ( fT9] l. ■ 
We note that the value in (fT9] l is approximately equal to 
0.25 when fc and p are large. 

Theorem |4] provides a hard guarantee on the worst-case 
system throughput; no matter what the delay offsets are, the 
system throughput is always larger than the value in ( fT9l l. 
Theorem |4] also indicates a tradeoff between the the lower 
bound and the sequence period. If we increase the value of 
fc, the sequence period is increased, but the lower bound on 
system throughput is also increased. 

We remark that the lower bound in Theorem |4] is not 
tight. The actual system throughput is higher than ( fT9] l. The 
next example compares the lower bound with the average 
throughput over relative delay offsets. 

Example 2. We consider an example with M 19 users, 
using CRT sequences with p = 37 and q = kp — 1. 
The throughput is plotted against the sequence period, while 
keeping the fraction of ones in each sequence fixed at 1 /p. This 
means that the fraction of time in which a user is transmitting, 
and hence the power of each user, is kept constant. We 
compare the lower bound in fT% with the average throughput 
obtained by simulation in Fig. [T] For each fc, 20000 delay 
offset combinations are randomly generated. The mean system 
throughput is about 0.31 In addition to the mean throughput, 
the maximum and minimum throughput obtained among these 
20000 delay offset combinations are also plotted. The variation 
of throughput diminishes as sequence period increases. The 
value of the lower bound ( fT9] l also increases. We see that the 
minimal observed system throughput is much higher than the 
lower bound in Theorem |4] 

For the shift-invariant protocol sequence set for nineteen 
users, the sequence period is 19^^, which is astronomical. 
Nevertheless, it has higher system throughput e^^ = 0.3679. 
For wobbling sequences, a lower bound of 0.25 system 
throughput [10, (6.12)] is guaranteed when the sequence 
period is 19* « 1.3 x 10^. From Fig. [T] a lower bound of 0.25 
can be obtained by using CRT sequences when the sequence 
period is about 1.1 x 10*, a roughly ten-fold reduction in 
sequence period. 




Max observed throughput 
Mean throughput 
Min observed throughput 
Lower bound 



Sequence Peirod 



Fig. 1. System Througliput of CRT Sequences for M = 19 users, p = 37. 



V. Conclusion 

In order to minimize the waiting time until a successful 
packet is sent, while maintaining a high level of system 
throughput, a class of protocol sequences with short period are 
constructed. After analyzing the crosscorrelation properties, 
we derive a lower bound on the system throughput. The con- 
structed sequences provides flexibility on the tradeoff between 
sequence period and the worst-case system throughput. 
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